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Abstract 

We consider nuclei produced in core collapse supernovae and subjected to a 
high neutron flux. We show that an accelerator mass spectrometry experiment 
that searched for traces of anomalous iron isotopes could set limits on the 
order of 10~ 20 — 10“ 25 on (or perhaps discover) Pauli principle violation by 
neutrons. A similar search for anomalous Co isotopes could set limits in the 
range 10~ 13 — 10~ 18 on Pauli principle violation by protons. We show that 
existing data on Oxygen can be used to set a limit of about 10 -1 ' in one 
proposed model of such violation. 


The Pauli exclusion principle, one of the the fundamental principles on which 
our basic understanding of the inner workings of systems ranging from atoms and 
nuclei to solids and liquids all the way to the stars and the universe rests, is clearly 
valid to a high degree of precision. Two simple examples in the domain of chemistry 
and physics are the patterns observed in the periodic table and the success of the 
shell model for nuclei. Until recently, however, there were hardly any quantitative 
tests of this fundamental principle. The best limit on possible deviation from the 
Pauli principle before 1987 was from an early experiment of Goldhaber and Scharff- 
Goldhaber [jj] in which beta rays from C 14 were allowed to impinge on lead and a 
search was made for K-shell x-rays. This bound on possible violation of the Pauli 
principle (characterized by a parameter f3 2 in this article and earlier ||) by electrons 
was at the level of 3%. In 1968, Fishbach, Kirsten and Schaffer 0 conducted a 
search for anomalous 9 Be and put an upper limit on the atmospheric density of 
such elements. Their result could be interpreted 0 to derive an upper limit of 
P 2 < 2 x 10" 3 . 

In 1987-90, several theoretical papers postulated ways to incorporate small 
deviations from Bose and Fermi statistics into quantum mechanics and held theories 
BI i 1- At this time, it was noted [§ that atomic spectroscopy could improve 
this bound for electrons to the level of 10 - '. Inspired by the theoretical interest in 
the subject, several dedicated experiments were planned around 1987-88 and have 
now considerably improved the limits for electrons B i I- These limits are now 
at the level of 1CU 26 - 10 2 '. As far as other particles go, strong limits for protons 


(though much weaker than for electrons) have been extracted by Plaga [|T(J from 
considerations of energy generation in the Sun. We have listed these limits in the first 
four rows of Table 1. Note that to date no limit has been established for neutrons. 
There also exist limits on deviations from bose statistics for bosonic systems such 


1 





as pions photons JTT| and spin zero atoms |T2[, which we do not address here. 

Meanwhile, on the theoretical front it was shown that deviation from Fermi 
and Bose statistics can be embodied in a satisfactory manner [6|] using the following 
“quonic” commutation relations p[] among the annihilation and creation operators 


(cp, a]): 


ctjoj — gajcq = Si 


( 1 ) 


The Hilbert space of states for this model is positive definite, though it leads to a 
non-local field theory |3] and a number operator which is non-polynomial. Also we 
note that for the fermionic particles, q = — 1 + /3 2 . 

In view of the fundamental role of the Pauli principle, it appears to us of 
considerable interest to expand the domain of our knowledge of the limits on Pauli 
violation (PV for short) to as many particles as possible. 

The purpose of this note is to point out potential, high precision tests of the 
Pauli principle, for neutrons and protons, making use of the high neutron and proton 
fluxes in a type II supernova and also to derive a bound from already existing limits 


on anomalous oxygen isotope by Hcmmick et al. At the risk of removing all 

suspense, we have included our results from this paper in Table 1. 
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Table 1. Pauli Principle Tests 


Particle 

Reference 

Method 

Upper limit on /3 2 

Electron 

Goldhaber and Golclhaber |I| 

K-capture 

0.03 

Electron 

Fishbach et al 0 

Anomalous 9 Be 

2 x 10" 3 

Electron 

Greenberg and Mohapatra f| 

Spectroscopy 

10~ 7 

Electron 

Kclleher 0 

Spectroscopy 

io - 7 

Electron 

Ramberg and Snow || 

Anomalous X-rays 

10~ 27 

Electrons 

Pomansky et al. 0 

Mass spectrometry 

o 

i 

to 

05 

Proton 

Plaga 0] 

Solar p-p cycle 

10- 15 

Neutron 

This paper 

Supernova Oxygen & 
Mass spectrometry 

10- 17 

(mo del- dep endent) 

Neutron 

This paper 

Supernova Iron & 
Mass spectrometry 

10~ 25 
is possible 

Proton 

This paper 

Supernova Cobalt & 
Mass spectrometry 

10~ 18 
is possible 


Our argument is as follows. About half of the terrestrial iron is produced during 
core collapse of massive stars to form Type II (or Type Ib/c) supernovae |L4|. There, 
nickel is produced and there is also a significant flux of free neutrons. After the 
supernova bounce, the ensuing shock wave ejects the nickel which decays to cobalt 
which, in turn, decays to iron. In the supernova explosion there is approximately a 
fraction of a second during which the nickel is formed and exposed to a high flux of 
neutrons with thermal energies of about 0.5 MeV. The exact value of the neutron flux 
depends on the electron fraction of the material ejected Y e , which in turn depends 
on the uncertain details of the supernova process. The other half of the terrestrial 
iron is produced in Type la supernovae by the same process, but the density is lower 
and the free neutron density is significantly lower since the supernova occurs as a 
thermonuclear explosion of a C+O white dwarf and there is little time for weak 
interactions to occur. Thus we will concentrate on iron produced in core collapse 
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supernovae. During the explosion, Pauli violation (PV) could permit the process 

n + A_1 Ni Ni* + 7 , (2) 

where A Ni* contains an extra neutron in an already occupied shell model state. For 
example it could have 3 neutrons in the 1 — s shell. Assuming that this (inner) 
configuration was not disturbed in the electron capture transitions from Ni* to Fe* 
(and did not inhibit those transitions), we would have today about the same relative 
abundance, up to a factor of two, of anomalous iron. Normal stable iron isotopes 
have A= 54, 56, 58. The unstable isotopes with A=55 and 59 are the products of 
shortlived isotopes and also have very short halflives. Hence they are not expected 
to occur in terrestrial radioactivities. Thus experimental detection of these two, 
which is well within the capabilities of accelerator mass spectrometry, would be an 
indication of Pauli violation. The abundances of 54 Fe and 58 Fe, from which 55 Fe* 
and 59 Fe* would be made, effectively, are 5.9% and 0.33%. 

We turn now to the question of what abundance ratio, Fe* /Fe, we can expect 
for a given amount of Pauli violation (PV). The probability that a neutron will be 
captured into a PV state is given by 

Py = /3 2 anvt (3) 


where f3 2 measures the suppression of the cross section a that would obtain in the 
absence of the Pauli principle, n is the neutron density, v is the velocity of the 
neutrons, and t is the time available. 

There should be three contributions to the cross section for the process of 
Eq(^[): magnetic dipole radiation from neutron spin flip; magnetic dipole radiation 
from Ni spin flip; and electric dipole radiation from Ni. Blatt and Weisskopf |15| 
give a useful discussion of capture and disintegration cross sections which we follow. 
We keep only the first contribution (neutron spin flip) since the others are smaller, 



and we attempt to be conservative in estimating the rate. We use the zero range 
approximation (Eq (4.27) of El)- We ignore possible scattering enhancement fac¬ 
tors, i.e. (1 — ay ) 2 where a is the scattering length and 7 = ( 2 pB) 1 ^ 2 with fi the 
n-Ni reduced mass, about the nucleon mass, and B the anomalous binding energy, 
on the order of 20 MeV. The capture cross section is then 

a = Tr(e 2 /hc)(hc/M n c 2 ) 2 (B/M n c 2 ){2B/E n ) 1/2 (4) 

= 10 " 3 O cm 2 , 

where E n = 0.5 MeV is the neutron energy. For the density of neutrons, we have 
attempted to be conservative in our estimate, and are guided by the results from 
assuming nuclear statistical equilibrium (NSE) (see Table 2 for some typical values 
in the region of the core that is ejected by the shock wave based on the model of 
Cooperstein, Kahana and Baron ||T 6 ||). 

I 11 the supernovae the material is not ejected with NSE abundances, but rather 
undergoes an alpha rich freeze out that has been suggested to be in quasi-statistical 
equilibrium |T7| (more alpha particles at a given (p,T,Y e ) than would be predicted 
by NSE). The neutron and proton densities drop precipitiously with temperature 
as the material expands. While the neutron flux is not determined by NSE, it is 
indicative. As the temperature drops during expansion, NSE is an increasingly poor 
assumption and so one has to both use a QSE code and know the expansion, T, and 
Ye history of the material, which is totally uncertain. Therefore the number X n used 
by us is simply a guess “guided” by the results in Table 2 which is based on assuming 
a reasonable equation of state for supernova material. We take p ~ 10 9 g cm -3 , and 
a neutron mass fraction X n ~ 10” 4 , giving n ~ 10 29 cm -3 . The thermal velocity of 
0.5 MeV neutrons is about 10 9 cm s _1 and the Ni is exposed to neutrons for about 
a 0.1 second. This gives 

P v ~ lO 8 /? 2 . 
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( 5 ) 


That is, we have 


P 2 < 10~ 8 EB 


( 6 ) 


where EB stands for an experimental upper bound on the relative abundance of 
anomalous iron. The above estimate is for the dominant isotope of iron i.e. 56 Fe. 
Since we are interested in 59 Fe*, we must multiply this probability by the factor 
3 x 10 -3 to take into account the low abundance of 58 Fe (which presumably implies 
the corresponding low abundance of the parent Ni isotope which participated in the 

, we 


PV reaction). Based on experience with such bounds for other elements [13 
can expect that limits in the range 10 ~ 13 to 10 -18 are readily achievable. Such an 
experiment would then set a PV limit for /3 2 on the order of 3 x 10 -19 to 3 x 10 -24 . 
If we use a similar argument for anomalous 55 Fe*, we would consider a parent Fe 
nucleus to be 54 Fe whose abundance is 5.9%. Thus, looking for 55 Fe* would provide 
bound of 1.5 x ICE 20 to 1.5 x ICE 25 . 


Table 2 . Illustrative Abundances for Supernova Conditions 


p (g cm 3 ) 

T (MeV) 

Y e 

X n 


v Q 

v Fe 

v Ni 

10 10 

1.0 

0.49 

5.2 (-3) 

1.8 (- 2 ) 

5.2 (-1) 

1,8 (- 1 ) 

5.0 (-4) 

10 10 

1.0 

0.495 

4.4 (-3) 

2.1 (- 2 ) 

5.5 (-1) 

1.6 (- 1 ) 

9.5 (-4) 

10 10 

1.0 

0.50 

3.7 (-3) 

2.5 (-2) 

5.7 (-1) 

1.4 (-1) 

1.6 (-3) 

10 9 

1.0 

0.49 

7.2 (-2) 

5.2 (-2) 

8.8 (- 1 ) 

5.6 (-10) 

7.3 (-14) 

10 9 

1.0 

0.495 

6.6 (- 2 ) 

5.6 (-2) 

8.8 (- 1 ) 

4.2 (-10) 

7.4 (-14) 

10 9 

0.5 

0.50 

1.1 (- 8 ) 

4.2 (-3) 

5.4 (-3) 

8.0 (-5) 

7.7 (-1) 

10 8 

1.0 

0.49 

3.0 (-1) 

2.8 (- 1 ) 

4.3 (-1) 

1.4 (-27) 

3.0 (-31) 

10 8 

1.0 

0.495 

2.9 (- 1 ) 

2.8 (- 1 ) 

4.3 (-1) 

1.3 (-27) 

3.1 (-31) 

10 8 

0.5 

0.49 

5.2 (-7) 

7.3 (-3) 

4.1 (-2) 

1.3 (-2) 

1.6 (- 1 ) 

10 8 

0.5 

0.495 

3.5 (-7) 

1.1 (- 2 ) 

4.2 (-2) 

4.8 (-3) 

2.8 (- 1 ) 

10 8 

0.3 

0.50 

0.0 

5.5 (-5) 

5.3 (- 6 ) 

8.7 (-10) 

9.9 (-1) 


We need to emphasize the fact that the estimate of neutron density in the 
relevant region is subject to considerable uncertainty. There is as yet no successful 
model of supernova explosion. The assumption of NSE is approximate, but should 
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not be unreasonable given that silicon burning to Ni is explosive and takes place at 
temperatures T > 0.5 MeV. In reality we would need to consider all (n, 7 ) reactions 
that occur in the nucleosynthetic process and we have tried to be conservative in our 
assumptions in order to derive a meaningful limit. The time available in the neutron 
bath and the density trajectory of the ejected material is a complicated function of 
explosion details. We believe however that these effects should not change the result 
quoted by more than a factor of 100. It is important to note that, even for the worst 
case scenario, we can expect a meaningful bound on (3 2 for neutrons from these 
considerations. 

Table 1 compares this result with other PV limits. The closest comparison 
is to the work of Plaga [[Hj who considered the effect of a possible triplet s-wave 
(hence PV) proton-proton bound state on the solar p-p cycle, achieving a bound of 
10~ 15 from the (assumed) approximate validity of the standard solar model. 

I 11 addition to the large neutron flux, the supernova has a large proton flux. 
Indeed, as noted above the proton flux is sufficiently large to overcome the coulomb 
barrier and convert most of the iron to nickel. This tells us that, for protons, the 
coefficient of (3 2 in Equation (3) must be on the order of one (or larger). As in the 
neutron case, we are assuming that the extra 1 -s proton neither inhibits the decays 
nor is disturbed by them. For the PV case, because the proton will be much more 
deeply bound, we might expect it to be considerably larger. Here we just assume 
it is of order 1. We thus see immediately that the bound on (3 2 is just equal to 
the bound on the abundance of the anomalous isotope, Co*. One could argue a 
10 3 enhancement in view of the fact that the abundance of Co* will be a fraction 
of the abundance of Fe while the abundance of normal Co is, at least on Earth, 
ICE 3 that of Fe. There are no long lived Co isotopes besides the stable 59 Co. The 
other isotopes are products of shortlived isotopes and also have very short halflives. 
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Moreover the decay energies of Co isotopes are much less than the 20 MeV or more 
that should be gained by putting the final proton in the Is state of Cobalt. Thus we 
expect that any detection of Z=27 nucleus with A 0 59 in a non-irradiated target 
would indicate Co* and hence Pauli violation. 

Finally, we note that existing data can be used to set a limit of interest for 
one PV model. One of us (RNM) has shown that if Pauli violation is characterized 
by the commutation relations (see ||) in Eq. (1), then a matrix representation for 
the annihilation operators is given by 


/ 0 0 0 ... \ 

10 0 .. 

0 0 0 .. 


(7) 


\. ... j 

with (3^ = 1 + q + ... + q n . For q close to —1, one sees that for odd n, /3 n is 
of the order 0 n+1 ^ 2 while, for even n, /3 n is of order 1. Thus, in a neutron rich 
environment, the second neutron is added relatively easily. The presupernova star 
has an oxygen shell with a neutron abundance around n rsj 10 20 . Eq (|) gives a 
capture probability of about 0.1 (3 2 . Middleton et al, some time ago, put a mass 
spectrometer limit on on anomalous oxygen, down to 20 AMU of about O* jO < 
1CT 18 which gives a PV bound, within the model of ||, of f3f < 10 _1 . We note 
however that the above matrix representation was derived in 0 for the case of a 
single oscillator. It is therefore not clear whether it will remain valid in an actual 
field theory generalization. One possible argument which may suggest its validity in 
a general field theory, can be given as follows. If we put the quon fields in a thermal 
bath, in the low temperature limit, we expect the partition function for the system 
to be dominated by the lowest energy level which will then be described by a single 




oscillator. For this case then the single oscillator representation should remain valid. 

In summary, we suggest use of mass spectrometry to limit Pauli principle vio¬ 
lations for neutrons and protons by exploiting the high nucleon fluxes in supernovae. 
A neutron limit of interest on a specific model is derived from existing data on Oxy¬ 
gen isotopes and much lower, model-independent limits from new experiments with 
iron and cobalt are possible. 
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